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ARTICLE INFO ABSTRACT 


Keywords: Swine acute diarrhea syndrome coronavirus (SADS-CoV) is a newly identified enteric alphacoronavirus that 
SADS-CoV causes fatal diarrhea in newborn piglets in China. Here, we propagated a virulent strain SADS-CoV/CN/GDWT/ 
Attenuation 


2017 in Vero cells for up to 83 passages. Four strains of SADS-CoV/GDWT-P7, -P18, -P48 and -P83 were isolated 
and characterized. Sequence alignments showed that these four novel strains exhibited 16 nucleotide mutations 
and resultant 10 amino acid substitutions in open reading frame 1a/1b, spike, NS3a, envelope, membrane and 
nucleocapsid proteins. Furthermore, a 58-bp deletion in NS7a/7b was found in P48 and P83 strains, which led to 
the loss of NS7b and 38 amino acid changes of NS7a. Pig infection studies showed that the P7 strain caused 
typical watery diarrhea, while the P83 strain induced none-to-mild, delayed and transient diarrhea. This is the 
first report on cell adaption of a virulent SADS-CoV strain, which gives information on the potential virulence 
determinants of SADS-CoV. 
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1. Introduction 


Swine acute diarrhea syndrome coronavirus (SADS-CoV), also 
known as porcine enteric alphacoronavirus (PEAV) or swine enteric 
alphacoronavirus (SeACoV), is a newly emerging porcine coronavirus 
and is the etiological agent of swine acute diarrhea syndrome (SADS), 
which can cause severe and acute diarrhea and high mortalities among 
neonatal piglets in China (Gong et al., 2017; Pan et al., 2017; Zhou 
et al., 2018a). SADS-CoV is an enveloped, positive and single-stranded 
sense RNA virus that belongs to genus Alphacoronavirus in the family 
Coronaviridae of the order Nidovirales (Zhou et al., 2018a). The SADS- 
CoV genome is approximately 27 kb and contains nine open reading 
frame (ORFs) that are ORFla, ORF1b, spike (S), envelope (E), mem- 
brane (M), nucleocapsid (N) and three accessory genes, NS3a, NS7a and 
NS7b (Zhou et al., 2018a; Wang et al., 2019). 

In 2017, SADS-CoV infection broke out in Guangdong Province of 
China and resulted in the death of around 25, 000 piglets within a few 
months (Zhou et al., 2018a). At that time, some farms used the virus 
solution collected from intestines of infected piglets to feed pregnant 
sows orally, to a certain degree which prevented the spread of SADS- 
CoV. However, SADS-CoV re-emerged in Guangdong on February 2019 
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and led to the death of around 2000 piglets in one farm (Zhou et al., 
2019a). So, there is an urgent need to isolate attenuated strains and 
develop vaccines to control SADS-CoV in pigs in China. Actually, the 
best way to prevent virus infection is vaccination of pregnant sows, 
which can induce lactogenic immunity, and their transmission to 
suckling piglets via colostrum and milk, are critical for early passive 
protection (Chattha et al., 2015). Like SADS-CoV, porcine epidemic 
diarrhea virus (PEDV) is also an enteric alphacoronavirus in pigs, but it 
can cause diarrhea to pigs of all ages and widespread in the world. To 
control PEDV, many live attenuated vaccines have been manufactured 
and licensed, for example, 83P-5 (Sato et al., 2011), attenuated CV77 
(Sun et al., 2016), DR13 (Park et al., 2012), and SM98 (Kim et al., 
2015). And more PEDV vaccine candidates including PEDVPT-P96 
(Chang et al., 2017), PC22A-derived strains (Lin et al., 2017), PEDV YN- 
144 (Chen et al., 2015), Zhejinag08 (Li et al., 2017), iPEDVPT-P96 (Kao 
et al., 2018), and VSVMT-SA19 (Ke et al., 2019) also have been de- 
veloped by serial cell passaging, cDNA clone construction or re- 
combinant attenuated virus vector. As a newly discovered alphacor- 
onavirus, recent studies regarding SADS-CoV have been mainly focused 
on molecular epidemiology, diagnosis, and experimental infection 
models (Fu et al., 2018; Zhoub et al. 2018, 2019; Xu et al., 2019), while 
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information on SADS-CoV attenuation and vaccine development is yet 
to be obtained. The emergence and re-emergence of SADS-CoV have 
posted a great threat to the pig industry in China, which make it ne- 
cessary to generate vaccines to prevent and control SASD-CoV infection. 

In the present study, we reported the attenuation of a virulent strain 
SADS-CoV/CN/GDWT/2017 via serial propagation in Vero cells for up 
to 83 passages. The pathogenicity of SADS-CoV/GDWT-P7 and -P83 
was evaluated in 6-day-old piglets and genetic variations related to 
attenuation were also identified by comparing full-genome sequences of 
SADS-CoV/GDWT-P7, -P18, -P48 and -P83. Our results indicated the 
virulent strain SADS-CoV/CN/GDWT/2017 was successfully attenuated 
and the SADS-CoV/GDWT-P83 strain could be further studied as atte- 
nuated vaccine candidate against SADS-CoV in China. 


2. Materials and methods 
2.1. Cells, animals and virus strain 


African green monkey kidney cells (Vero cells) were kept in our 
laboratory and cultured in Dulbecco's modified Eagle's medium 
(DMEM, HyClone) supplemented with 10% fetal bovine serum (FBS, 
HyClone) at 37°C within 5% CO, incubator. The SADS-CoV/CN/ 
GDWT/2017 strain (GenBank accession number MG557844) was iso- 
lated from diarrhea samples of suckling piglets in a commercial pig 
farm in Guangdong Province in April 2017 (Zhou* et al., 2018) and was 
kept in our laboratory. 26 five-day-old healthy piglets that didn't suck 
colostrum after birth were purchased form a commercial farm in 
Guangdong with no history of SADS. All piglets were tested negative for 
SADS-CoV, PEDV, porcine deltacoronavirus (PDCoV), pseudorabies 
virus (PRV), rotavirus (RV) and transmissible gastroenteritis virus 
(TGEV) by virus-specific qPCRs on rectal swabs, and were detected 
antibody negative for porcine reproductive and respiratory syndrome 
virus (PRRSV) (Li et al., 2013; Mai et al., 2017; Wu et al., 2017; Zhou 
et al., 2019a,b). 


2.2. Serial passage of the virus 


Vero cells were cultured to 80%—-90% confluence and washed three 
times with phosphate-buffered saline (PBS) (20% w/v). Then, cells 
were infected with the SADS-CoV/CN/GDWT/2017 strain and in- 
cubated with serum-free DMEM containing 10 ug/ml trypsin (Gibco). 
Infected cells were maintained at 37 °C with 5% CO. and monitored 
daily for cytopathic effects (CPE). When 90% of the cells showed CPE, 
cells were harvested, frozen and thawed three times, and centrifuged 
for 10minat 8000 g. The supernatant was collected and stored at 
— 80°C as the stock for the next passage. Using the same manner, 82 
subsequent passages were conducted in Vero cells. Beginning at the 8th 
passage in cell culture, the virus was plaque purified every ten passages. 
The viral stocks every five or ten generations were titrated by per- 
forming a 10-fold serial dilution of supernatants in 96-well plates in 
triplicate per dilution, and the virus titration was determined by TCIDso 
method (Reed and Muench, 1938). 


2.3. Genome sequencing, alignment and phylogenetic analysis 


The complete genomes of SADS-CoV isolates were amplified using 
previously reported primer pairs (Zhou et al., 2019a,b). PCR assays 
were performed with the following thermal profile: 95 °C for 3 min, 30 
cycles of 95 °C for 30s, 55°C for 30s, and 72°C for 60s, followed by a 
final 10 min extension at 72°C. The products were purified following 
the manufacturer's instructions of Gel Band Purification Kit (Takara, 
China) and then cloned into the pMD19-T vector (Takara, China) and 
transformed E. coli DH5a competent cells. The positive clones were 
screened out and sent to Beijing Genomics Institute (Shenzhen, China) 
for further sequencing. The nucleotide sequences were assembled and 
aligned using the DNASTAR program (DNAStar V7.1, Madison, WI, 
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USA). Phylogenetic trees were constructed using the neighbor-joining 
method in MEGA 7.0 software with bootstrap analysis of 1000 re- 
plicates. Percentages of replicate trees in which the associated taxa 
clustered are shown as nearby branches (Chenna et al., 2003; Tamura 
et al., 2004; Kumar et al., 2016). 


2.4. Animal experiment design 


26 piglets were randomly divided into three groups including two 
SADS-CoV- infected groups each containing nine piglets and one control 
group with eight piglets. After one day of acclimation, piglets in two 
inoculated groups were orally with 12 mL of 10° TCIDs9/mL of SADS- 
CoV/CN/WT/2017-P7 and -P83, respectively. Piglets in the control 
group were mock-inoculated in parallel with 12mL of DMEM. Fecal 
consistency and body weight for each piglet were recorded or scored 
daily. Rectal swabs were also collected every day to detect viral shed- 
ding. When diarrhea was observed in infected groups, three piglets in 
each group were randomly selected to euthanized and necropsied to 
collect tissue samples for further analyses. 


2.5. RNA extraction and qRT-PCR analysis 


SADS-CoV RNA was extracted from fecal or tissue samples following 
the manufacturer's recommendations of AxyPrep’™ Body Fluid Viral 
DNA/RNA Miniprep Kit (Axygen Scientific, Inc). Reverse transcription 
was carried out to synthetized cDNA according to the manufacturer's 
protocol of Prime ScriptTM RT Reagent Kit with gDNA Eraser (Takara, 
China). The TaqMan-based qRT-PCR assay established by our lab (Zhou 
et al., 2018b) was employed to amplify cDNA with specific primers 
(forward: 5’ - CAGGTCTTGGTGTTCGCAATCG-3’; reverse: 5’- ACCGTG 
CTGAACGAGGTCACT-3’), and probe (5’-FAM- TCACAGTCTCGTTCTC 
GCAATCA-TAMRA-3’) under the following protocol: 95°C for 30s for 
initial denaturation followed by 45 cycles of 95 °C for 5s and 62°C for 
30s. Ten-fold serial dilutions of the pMD19-T-SADS-CoV-N plasmid 
were used to generate a standard curve. The viral genomic copy num- 
bers were calculated by comparing 1ul of cDNA with a standard 
sample. 


2.6. Histopathology and immunohistochemical staining 


Tissue samples including duodenum, jejunum and ileum collected 
from the challenged and mock-infected piglets were fixed in 10% for- 
malin for 24 h at room temperature, then dehydrated in graded ethanol, 
embedded in paraffin, and cut to 3 um sections (Leica, Germany). The 
paraffin-embedded samples were deparaffinized in xylene, and washed 
in decreasing concentrations of ethanol. After that, the deparaffinized 
intestinal tissue sections were stained with hematoxylin and eosin 
(Sangon, China) for histopathology or subjected to im- 
munohistochemistry assay using SADS-CoV-specific mouse antisera 
(1:400) and goat anti-mouse IgG (H + L)/HRP (1:50) through proce- 
dures described previously with some modifications (Zhou et al., 
2018a). 


3. Results 


3.1. Biological and growth characteristics of the virus strain in serial of 
passages 


During the serial propagation, viral titers of passage (P) 12, 23, 33, 
38, 48, 53, 58, 63, 68, 73 and 83 were measured. The result showed 
that infectious titer of the passaged virus ranged from 10* to 107° 
TCIDs9/mL and rose sharply before P38, then slowed down and reached 
the highest titer at P83 (Fig. 1). As the same time, with the increase of 
the passage, CPE appeared more quickly. In P7-infected Vero cells, CPE 
appeared on 24h post infection (hpi) including syncytial and vacuole 
formation, while over 90% of P83-infected cells died and floated at this 
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Fig. 1. Viral titers of SADS-CoV/CN/GDWT/2017 strain propagated in Vero 
cells after serial passage. 


timepoint (Fig. 2). These results indicated an increasing adaptability for 
the virus in cell serial passage. 


3.2. Sequence analysis 


Four new complete genomes of SADS-CoV/GDWT-P7, -P18, -P48 
and -P83 strains (accession number MK994934 to MK994937) were 
sequenced and available in GenBank. The length for the complete 
genomes of SADS-CoV/GDWT-P7 and -P18 stains was 27165 bp, and 
that of SADS-CoV/GDWT-P48 and -P83 was 27107 bp. Sequence 
alignments showed nucleotide identities between all available SADS- 
CoV genomes ranged from 99.3% to 100%. The four novel full-length 
sequences shared high nucleotide identities of 99.9%-100% with each 
other and 99.5%-100% with other reported complete genomes of 
SADS-CoV. Phylogenetic analyses based on complete genomes demon- 
strated that all SADS-CoV strains clustered together to form a well-de- 
fined group and separated from other alphacoronaviruses (Fig. 3). 

After comparing full-length sequences of four new strains with that 
of SADS-CoV/CN/GDWT/2017, no insertion events were observed, but 
16 nucleotide (nt) point mutations and the resultant 10 amino acid (aa) 
substitutions were found (Table 1). The ORFla/1b had most mutations 
containing 9nt and 5 aa changes, followed by S and N gene each har- 
boring 2nt and 1 aa changes. The E, M and NS3a had one aa sub- 
stitution, respectively. Besides these mutations, a large deletion of 58 
bp was found in SADS-CoV/GDWT-P48 and -P83. This 58 bp-deletion 
led to the loss of initial codons for both NS7a and NS7b and a reading 
frame shift for NS7a. As a result, the NS7b gene was lost in genomes of 
these two high-passage strains. For the Ns7a gene, the reading frame 
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moved forward and caused an overlap of 107bp between the N and the 
NS7a. Finally, 20 additional aa were encoded and 18 aa substitutions 
were substituted in the initial part of NS7a (Fig. 4). 


3.3. Pathogenicity difference between the low- and the high-passage strains 


To explore whether these mutations and deletions of amino acids in 
high-passage strains can lead to virulence changes, SADS-CoV/GDWT- 
P7 (group A) and SADS-CoV/GDWT-P83 (group B) at same TCIDs9 
value were inoculated into six-day-old piglets. Clinical scoring of fecal 
consistency in challenged and control groups were shown in Table 2. At 
one day post-inoculation (DPD, all piglets (9/9) in the SADS-CoV/ 
GDWT-P7-inoculated group presented with typical SADS-CoV-asso- 
ciated, loose-to-watery diarrhea. And the diarrhea peaked on 2-3 DPI 
(score => 2), accompanying with watery feces ejecting from wet, red and 
swollen anus. Piglets were dispirited, not moving and even vomiting. 
The clinical signs lasted for six days. At 7 DPI, piglets began to recover 
and diarrhea disappeared. On the contrary, only two piglets (B3 and 
B6) in the SADS-CoV/GDWT-P83-inoculated group developed loose-to- 
semi-fluid diarrhea (score<2) on 5 DPI, and rapidly recovered on 7 
DPI. The other piglets in this group showed no obvious clinical symp- 
toms. In the mock-infected group, no clinical signs of diarrhea were 
observed. 

The mean body weight in each group was investigated to assess 
clinical effects of diarrhea resulted by SADS-CoV/GDWT-P7 and -P83 
infections in piglets. As shown in Fig. 5, piglets in the control group 
obviously weighed more than piglets in two challenged groups. Con- 
trary to mock-infected piglets continuously gaining the body weight, 
piglets in SADS-CoV/GDWT-P7-infected group almost stopped growing 
on 1-2 DPI. After that the body weight began to slowly increase. For 
piglets inoculated by the SADS-CoV/GDWT-P83 strain, the body weight 
was basically same with that of the control group within three days 
after challenging, but nearly didn't alter in the following days, and on 5 
DPI started to increase again. 

The results of rectal swabs collected from SADS-CoV/GDWT-P7-in- 
noculated group showed that SADS-CoV was detected positively from 2 
to 7 DPI, which was approximately consistent with that of diarrheal 
clinical signs. Although the moderate diarrheal symptom of a piglet 
(A9) disappeared after 3 DPI, SADS-CoV was still detected positively 
from the rectal swab of this piglet until 7 DPI. In the SADS-CoV/GDWT- 
P83-inoculated group, rectal swabs from six piglets were detected po- 
sitive for SADS-CoV, and the number was more than that of piglets with 
clinical signs of diarrhea. The piglet numbered B3 with the moderate 
diarrhea was detected positive for SADS-CoV from 2 to 5 DPI; it was the 
longest duration of viral shedding in this group. Three piglets including 
two euthanized were detected negative for SADS-CoV during the ex- 
periment. All piglets in the control group were also negative for SADS- 
CoV (Table 3). 


24h 48h 


Fig. 2. Cytopathic effects (CPE) caused by SADS-CoV/CN/GDWT/2017-P7 or -P83 (100 x ). Vero cells were seeded into T-25 flasks and infected at a multiplicity of 


infection (MOI) of 0.01. 
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Fig. 3. Phylogenetic analysis of the complete genome of SADS-CoV and reference coronavirus species. The tree was constructed using MEGA 7.0 software with 
neighbor-joining methods and 1000 replicate sets on bootstrap analysis. Four new complete genomes sequences studied in this work were indicated with “black 


blots”. 
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Table 1 
Changes in nucleotides and amino acids in SADS-CoV/GDWT-P7, -P18, -P48 and -P83 strains. 


Gene Nucleotide Change Amino Acid Substitution 
Position SASADS-CoV/CN/GDWT/2017 P7 P18 P48 P83 
OFR1la NSP2 683 C C T T T A17V 
NSP3 3838 G G G T I A131S 
NSP3~NSP4 6230 T T T C C L=S 
NSP4 7832 A A A G G D47G 
NSP6 10242 C C C T T - 
NSP6 10737 A A C C C - 
NSP6 10831 G G G A A V266I 
NSP9 11994 T T T G G - 
OFR1b NSP10~NSP11 16763 C C C T T - 
S S1 21166 T T T T C - 
S2 22497 A A A C C E671A 
NS3a 23987 C A A A A P37H 
E 24678 G G G T T R44M 
M 24832 C C TE di f L17F 
N 26131 C C C T T A217V 
26447 C C C C T - 


A EEA TODAAROOTEUAROOTOMAROGTO AAG NODOUAOAOTOUARNOOUEGOGY 


Majority ATTAATTAGTTCAACTAAACATGAACCAGCTTATCTTTTTCCTGTGTATGATGTGCTGTTATGCTATTTTTTTTGATTGGCT 

— se a ee ee eC OCS 

26610 26620 26630 26640 26650 26660 26670 26680 

ee ee ee ee ee ee ee ee 

SADS-CoV-GDWT-P7.seq ATTAATTAGTT CAACTAAACAT GAACCAGCTTATCTTTTTCCTGIGIATGATGIGCTGITATGCTATTTTTTTTGATTGGCT 

SADS-CoV-GDWT-P18.seq ATTAATTAGTTCAACTAAACATGAACCAGCTTATCTTTTTCCTGTGTATGATGTGCTGTTATGCTATTTTTTTTGATTGGCT 

SADS-CoV-GDWT-P48.seq RT irene eareesiiicn ens annie mamCEIeERNR EM Mone TTGGCT 

SADS-CoV-GDWT-P83.seq RAI iit sce ea cern aaa al ele TTGGCT 
B 


[<x] Consensus MLMPQCSLRHLQLLOF PKI LLNLLI RLLWAMRLLAI STDWLFNLFFYACQVNTWOEFAFSCNWSWGLFFEDFSTWEKCLSVVLI GTI AAASFMFADFAVEVFDLFERFFI NVGRFCRFV- 
4 Sequences 10 20 60 70 80 90 100 110 120 
7 NS7a.seq 
18 NS7a.seq 
48 NS7a.seq MLMPQCSLRHLQLLOFPKI LLFLRRRLRWAMRLLI SSTNWLFNLFFYACQVNTWOEFAFSCNWSWSLFFEDFSTWEKCLSVVLI GTI AAASFMFADFAVEVFDLFERFFI NVGRFCRFV. 
83 NS7a.seq MLMPQCSLRHLOQLLOFPKI LLFLRRRLRWAMRLLI SSTNWLFNLFFYACQVNTWOEFAFSCNWSWSLFFEDFSTWFEKCLSVVLI GTI AAASFMFADFAVEVFDLFERFFI NVGRFCRFV. 


Fig. 4. Alignments of NS7a/7b genes of SADS-CoV/GDWT-P7, -P18, -P48 and -P83. A, the green solid boxes indicate the location of a 58 bp-deletion in P48 and P83 
strains. B, the light blue regions indicate amino acid changes of NS7a proteins of P48 and P83 strains. (For interpretation of the references to colour in this figure 
legend, the reader is referred to the Web version of this article.) 


Table 2 
Clinical scoring of diarrhea. The scores were graded as follows: 0, normal; 1, loose diarrhea; 2, semi-fluid diarrhea; and 3, watery diarrhea. DPI, days post- 
inoculation. 


DPI SADS-CoV/GDWT-P7-inoculated group SADS-CoV/GDWT-P83-inoculated group Mock-infected gourp 
Al A2 A3 A4 A5 A6 A7 A8 A9 Bl B2 B3 B4 B5 B6 B7 B8 B9 Cl C2 C3 C4 C5 C6 C7 C8 

0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 O 0 oO O 0 oO O 
1 3 3 1 3 3 2 1 2 1 0 0 0 0 0 0 0 0 0 0 O 0 oO O 0 oO O 
2 3 3 2 3 3 3 2 3 2 0 0 0 0 0 0 0 0 0 0 O 0 oO O 0 oO O 
3 3 3 3 3 3 1 2 2 1 0 0 0 0 0 0 0 0 0 0 O 0 oO O 0 oO O 
4 3 / 2 2 / 0 1 ff 0 i 0 0 0 7 0 rs 0 0 0 / 0 oO / 0 oO / 
5 2 7 1 1 7 0 0 "4 0 / 0 2 0 7 1 vA 0 0 0 / 0 oO / 0 oO / 
6 1 7 0 0 7 0 0 jf 0 / 0 1 0 / 0 / 0 0 0 / 0 oO / 0 oO / 
7 0 i 0 0 i 0 0 jf 0 / 0 0 0 / 0 / 0 0 0 / 0 oO / 0 oO / 
8 0 / 0 0 / 0 0 iP 0 4 0 0 0 / 0 i 0 0 oO F 0 oO / 0 oO / 
9 0 7 0 0 7 0 0 7 0 / 0 0 0 / 0 / 0 0 0 / 0 oO / 0 oO / 


/, A2, A5, A8, B1, B5, B7, C2, C5 and C8 were euthanized and necropsied for collecting tissue samples. 
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Fig. 5. Daily change in the body weight of 6-day-old piglets after SADS-CoV/ 
GDWT-P7 or -P83 inoculation. The daily differences in mean body weight for 
each group were shown as mean + standard error of the mean (SEM). 


3.4. Tissue tropism 


Twelve tissue samples including heart, liver, spleen, lung, kidney, 
tonsil, stomach, duodenum, jejunum, ileum, mesenteric lymph node 
and inguinal lymph node were collected to determine the different 
tissue tropism between the low- and high-passage. For the low-passage, 
P7 strain, viral RNA could be detected in samples of tonsils, stomach, 
duodenum, jejunum, ileum and mesenteric lymph nodes, among which 
the ileum had the highest viral load of 4.7710g10 copies/uL. However, 
the viral RNA of the SADS-CoV/GDWT-P83, was only detected in ileum 
and mesenteric lymph nodes and with same viral copies of 1.67log10 
copies/uL (Fig. 6). 


3.5. Gross pathology, histopathology, and immunohistochemistry 


The intestinal wall of small intestine of piglets inoculated by the 
SADS-CoV/GDWT-P7 strain became thinner, permeable, and the in- 
testine was filled with a large number of yellow water feces. The di- 
gestive contents in the cecum and colon were also watery and without 
pigmentation. Mesenteric congestion was obvious and accompanied by 
individual bleeding points. While in SADS-CoV/GDWT-P83-inoculated 
piglets, thick intestinal wall, normal mesentery and pigmentation were 
observed, which was almost similar with those of mock-infected piglets 
(Fig. 7). 

Histological changes associated with SADS-CoV infection were ob- 
served in the duodenum, jejunum and ileum samples collected from 
SADS-CoV/GDWT-P7-inoculated group. As shown in Fig. 8, SADS-CoV/ 
GDWT-P7-infected piglets had clear histologic lesions with villus 
blunting and atrophying, nuclei dissolving and shedding, mucosal in- 
flammatory cells increasing and sever adhesion between villi. While the 
normal villous epithelium was overserved in piglets in SADS-CoV/ 
GDWT-P83-inoculated group and the mock-infected group. Consistent 
with the histological results, SADS-CoV antigen signals appeared brown 
and were detected in epithelial cells in the duodenum, jejunum and 
ileum of SADS-CoV/GDWT-P7-inoculated piglets. No SADS-CoV anti- 
gens were detected in samples of piglets inoculated with SADS-CoV/ 
GDWT-P83 or mock-inoculated piglets (Fig. 9). 


4. Discussion 


As a newly emerging porcine enteric coronavirus, research on SADS- 
CoV attenuation has not been reported yet. Considering its high 
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pathogenicity, there is an urgent need to isolate attenuated strains and 
develop efficient vaccines to control SADS-CoV. In the present study, 
for the first time we serially propagated the virulent strain SADS-CoV/ 
CN/GDWT/2017 in Vero cells for up to 83 passages. Four novel strains 
SADS-CoV/GDWT-P7, -P18, -P48 and -P83 were successfully isolated 
and their genetic and pathogenic characteristics were ascertained. 
Compared to SADS-CoV/GDWT-P7-inoucalted piglets, SADS-CoV/ 
GDWT-P83-inoculated piglets showed a mild and transient diarrhea 
with shorter viral shedding, weakened tissue tropism and undetected 
histopathological lesion, which provides evidence of viral attenuation 
and highlights SADS-CoV/GDWT-P83 strain as a potential, live-atte- 
nuated SADS-CoV vaccine candidate. 

Entire genomes of these four isolated strains were sequenced to 
decode genetic changes that emerged during serial cell passage. By 
comparing with the sequence of SADS-CoV/CN/GDWT/2017, nucleo- 
tide mutations in ORFla/b, N, S, E, M and NS3a were identified. Our 
result showed that half of the non-silent mutations (5/10) were located 
in the ORF1a gene, i.e. NSP2, NSP3, NSP4 and NSP6. It is well known 
that the NSPs of coronaviruses are associated with viral replications, 
many of which are considered to be interferon antagonists (Xing et al., 
2013; Chang et al., 2017). Previous findings of PEDV attenuation also 
revealed relatively large nucleotide changes in ORFla/1b (Chen et al., 
2015; Lee et al., 2017), which suggests this may be a common feature 
for both SADS-CoV and PEDV. As a major surface protein, the cor- 
onavirus S protein possesses two subunits S1 and S2, which are re- 
sponsible for cell receptor interaction and cell membrane fusion, re- 
spectively (Hulswit et al., 2016). In the present study, only one 
nucleotide mutation and resultant one amino acid substitution was 
found in the S2 region, and no insertion or deletion was observed in S 
gene, which were different with the discovery of PEDV. Some cell 
culture-attenuated PEDV strains, e.g., PEDV YN-144 (Chen et al., 2015), 
PEDVPT-P96 (Chang et al., 2017) and PC22A-derived strains (Lin et al., 
2017), had many or most of amino acid changes in the region of S gene. 
Insertion mutations that resulted in four amino acid changes were 
found in the S1 region of cell culture-passaged PEDV CH/HNPJ/2016 
strains (Liu et al., 2019). However, there was also no nucleotide change 
observed in S gene of attenuated PED-CUP-B2014 strain (Park et al., 
2018). These results indicated multiple mutation patterns of virus at- 
tenuation during cell passaging. For SADS-CoV, more analyses are 
needed to explore potential molecular mechanisms of attenuation to 
provide valuable data for vaccine developments. 

Besides structural and non-structural proteins, the genome of SADS- 
CoV also encodes three accessory genes, NS3a, NS7a and NS7b (Zhou 
et al., 2018a; Wang et al., 2019). NS3a locates between S and E protein, 
corresponding to ORF3 in PEDV genome or ORF3a/b in TGEV genome 
(Zuniga et al., 2016; Wang et al., 2019). Many groups have demon- 
strated that cell-adapted PEDV or TGEV strains display either internal 
truncation or amino acid changes in the ORF3 gene, which suggests that 
ORF3 is of key importance for virus adaption and virulence attenuation 
(Zhang et al., 2007; Park et al., 2012; Jengarn et al., 2015; Wongthida 
et al., 2017). Similar to these findings, we observed the amino variation 
in the NS3a region of SADS-CoV/GDWT-P7 and other three new iso- 
lated strains, and we also noticed that the pathogenicity of P7-infected 
piglets was lower than that of the wild type SADS-COV strain. Although 
presenting typical SADS-CoV-associated diarrhea, P7-innoculated pig- 
lets survived and recovered in a week. The amino variation of NS3a in 
our assay seems to contribute to the virulence attenuation, which need 
further analysis to confirm the link between NS3a and attenuation of 
SADS-CovV. 

During serial cell passage of the present study, the largest nucleotide 
changes, a 58 bp-deletion, occurred in NS7a and NS7b of SADS-CoV/ 
GDWT-P48 and -P83 strains. This deletion led to the loss of NSA7b 
protein expression and the 38 amino acid-change of NS7a protein. 
Generally, the deletion of accessory genes could result in attenuated 
viruses, such as ORF7 deletion in TGEV (Ortego et al., 2003) and 
ORF7a/b deletion in feline infectious peritonitis virus (FIPV) (Haijema 
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Fig. 6. Copy numbers of SADS-CoV genomes in different tissues of piglets in- 
fected by SADS-CoV/GDWT-P7 or -P83 strains at 4 DPI. 


et al., 2004); but there was also a report that the virulence of ORF7- 
deleted TGEV was enhanced compared with that of parental viruses 
(Cruz et al., 2011). For SADS-CoV, as the total deletion of 7b protein 
and large substitutions of 7a protein appeared in the middle and late 
stage of cell passage, as well as the P83 strain presenting low virulence 
to piglets, we speculated that these changes in NS7a/7b may also 
contribute to the attenuation of SADS-CoV. In addition to virulence 
modulation, protein 7 also plays a role in interactions between virus 
infection and innate immunity. The absence of TGEV protein 7 in- 
creased expressions of IFN-f, IFN-stimulated genes and_proin- 
flammatory genes (Cruz et al., 2013). FIPV 7a protein acted as a type I 
IFN antagonist, when cooperating with ORF3-encoded proteins 
(Dedeurwaerder et al., 2014). Functions of protein 7 of SADS-CoV is 
necessary to study further to clarify its roles in cell adaption and 
virulence alteration. Collectively, our results showed that all nine ORFs 
of SADS-CoV exhibited nucleotide and amino acid changes during serial 
propagation, which indicates the attenuation of virus may be multi- 
factorial. These variations mentioned above along the genome may 
constitute a combination to have an impact on cell adaptation and 
virulence attenuation of SADS-CoV, or it is possible that some in- 
dividual mutations may have key determinants, all of which need more 
investigations by the reverse genetic method. 
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5. Conclusions 


In summary, SADS-CoV/CN/GDWT/2017 was successfully atte- 
nuated through serial passaging in Vero cells as evidenced by enhanced 
virus titers in cells and reduced pathogenesis in neonatal piglets. 
Genetic investigations revealed that 10 non-silent mutations in ORFla/ 
1b, S, E, NS3a, M and N, and a 58-bp deletion in NS7a/7b contributed 
to cell adaption and virulence attenuation. To our knowledge, this is the 
first report to serial passage culture and attenuation investigation of a 
virulent SADS-CoV strain. The results of this study are important for 
mapping the virulence determinants of SADS-CoV and for the devel- 
opment of effective vaccines against SADS in China. 
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Fig. 7. Intestinal dissection of experimental piglets infected by SADS-CoV/GDWT-P7 or -P83 strains at 4 DPI. Thin-walled intestinal tract of containing yellow water 
feces was indicated by the arrow. (For interpretation of the references to colour in this figure legend, the reader is referred to the Web version of this article.) 
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Fig. 8. H&E-stained intestinal tissue section of experimental piglets infected by SADS-CoV/GDWT-P7 or -P83 strains at 4 DPI. Blunt intestinal villus was indicated by 
arrows. 
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Fig. 9. Immunohistochemistry (IHC) of the small intestines of piglets inoculated with SADS-CoV/GDWT-P7 or -P83 strains at 4 DPI. Arrows indicate the antigen 
signals of SADS-CoV. 
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Abbreviations 


SADS-CoV: swine acute diarrhea syndrome coronavirus 
SADS: swine acute diarrhea syndrome 

PEAV: porcine enteric alphacoronavirus 

SeACoV: swine enteric alphacoronavirus 

ORF: open reading frame 

S: spike 

E: envelope 

M: membrane 

N: nucleocapsid 

PEDV: porcine epidemic diarrhea virus 

PDCovV: porcine deltacoronavirus 

PRRSV: porcine reproductive and respiratory syndrome virus 
PRV: porcine pseudorabies virus 

RV: rotavirus 

TGEV: transmissible gastroenteritis virus 

PBS: phosphate-buffered saline; 

CPE: cytopathic effects 

hpi: h post infection 

DPI: day post-inoculation 


